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Abstract 

The surface stabilized MnO(100)-like monolayer, characterised by a regular c(4x2) distribution 
of Mn vacancies, is studied by hybrid functionals and discussed in the light of available scanning 
tunneling microscopy and high-resolution electron energy loss spectroscopy data. We show that the 
use of hybrid functionals is crucial to account for the intermingled nature of magnetic iner actions, 
electron localization, structural distortions and surface phonons. The proposed Pd(100) supported 
MJI3O4 structure is excellently compatible with the experiments previously reported in literature. 
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I. INTRODUCTION 



The growth of oxide ultrathin films on metal substrates represents a fascinating and 
active field of research in modern surface science, primarily because of the particular and 
novel chemical and physical properties, substantially different and to a certain extent richer 
than the corresponding bulk counterparts. The increasing interest stems from the numerous 
technological applications and the challenging fundamental understanding that have boosted 
intensive experimental and theoretical efforts. 

In this context, manganese oxides represent a distinguished example, partly because of 
the complex interplay between orbital, spin and lattice degrees of freedom which lead to 
prominent phenomena such as colossal magnetoresistence, metal-insulator transitions and 
exotic magnetic behaviours |l|, [2], [], 0, [5], but also from a chemical point of view considering 
the wide range of applications comprising catalysis, environmental waste treatments, pro- 
duction of water purifying agents and a^ne/dry-^ battertj fl Q . As such, thin 
films of manganese oxides on metal supports are promising candidates with potential new 
applications in many different areas. 

Thin films of MnO have been successfully prepared on various noble metal substrates, 
such as Ag(001)[10], Rh(OOl) 11] and Pt(lll) 12, 13]. In a series of previous studies we have 



ocusing on both 
In 



15, m 



provided a detailed description of the Mn oxide phases formed on Pd(100 
high [20-30 monolayers (ML)] and low (below 1 ML) coverage regimes [14 
particular, we have shown that upon deposition of 20-30 ML epitaxial MnO(100) films can 
be grown which can be preferentially converted either into MnO(lll) or Mn 3 O4(001) by ap- 
propriate tuning of temperature and oxygen pressure. Below 1 ML, a complex surface phase 
diagram was reported, where nine different novel Mn oxide phases have been detected which 
belong to two distinct oxygen pressure regimes and are characterised by well-defined struc- 
tural and vibrational properties, as established by scanning tunnelling microscopy (STM), 
low-energy electron diffraction (LEED), high- resolution electron energy loss spectroscopy 
(HREELS) and density functional theory (DFT). 

Here, by means of first principles calculations, we aim to explore the structure and prop- 
erties of one particular Mn x O y -Pd(100) interfacial phase, namely the c(4x2)-Mn3C>4 struc- 
ture which is formed at intermediate O2 pressures (around 5xl0 -7 mbar). Interestingly, 
this c(4x2)-Mn3O4/Pd(100) structure appears to be very similar to the recently observed 
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c(4x2)-Ni3C>4 phase obtained upon reactive evaporation of nickel on Pd(100), which has been 
identified as a compressed substoichiometric NiO(lOO) monolayer characterised by a regular 
rhombic distribution of Ni vacancies 0,0,0- The analogy between Pd(100)-supported 
MH3O4 and Ni 3 04 c(4x2) structures is likely to originate from the similarity between the 
corresponding bulk parent compounds. MnO and NiO crystallise in the same rhombohe- 
drally distorted fee structure, display a type II antiferromagnetic spin ordering, and the 
electronic nature is dominated by the partially filled localised 3d states. However, a major 
difference is observed in the lattice constants of MnO (aM n o = 3.14 A) and NiO (aNio=2.95 
A) that renders the analogy between M113O4 and M3O4 c(4x2)/Pd(100) hardly predictable. 
Specifically, for the deposition of MnO and NiO on Pd(100) (apd=2.75 A), this difference 
in the lattice constants results in a much higher positive lattice mismatch for MnO (14 %) 



than for NiO (7 %). Despite that, bulk-like MnO(lOO) pj and NiO(lOO) [2l[ films has 
been both successfully grown on a Pd(100) substrate. Summing up, in spite of the different 
lattice mismatches, the strong similarities of MnO and NiO seem to support the first glance 
analogy between the corresponding interfacial Pd(100) supported c(4x2) phases. 

In the present work we address this issue and interpret the Mn304 structure as a com- 
pressed epitaxial MnO(lOO) monolayer with a rhombic c(4x2) array of manganese vacancies. 
Our study is mainly concerned with the ab initio investigation of structural, electronic, vi- 
brational and magnetic properties, but a link to the experimental findings will be recalled 
when necessary. 

The paper is organised as follows: In Sec. [Ill we describe the computational tools, whereas 
the results will be presented and discussed in Sec. IIII1 Finally, in Sec. HV]we draw conclusions. 

II. COMPUTATIONAL DETAILS 

The results presented in this work were obtained using the Vienna ab initio simulation 



package (VASP) 22] within standard (Kohn Sham theory) and generalised (hybrid) density 



functional theory (DFT) 23|, |24|. All calculations have been performed using the projector- 
augmented-wave (PAW) method within the generalised gradient spin density approximation 



to the DFT in the Perdew-Burke-Ernzerhof parametrisation scheme 



satisfactory application of the HSE (Heyd-Scuseria-Ernzerhof) 



28, 2J] 



271 ] . Prompted by the 



hybrid DF' 



on the physical properties of manganese oxides presented in a series of recent papers 



scheme 

3,B,Q, 



.3 



we have adopted the same hybrid formalism in the present work to calculate the electronic 
properties and phonon frequencies of the most favourable models. 

Unlike standard DFT, HSE employs an admixture of Hartree-Fock (HF) and PBE ex- 
change in the construction of the many body exchange (x) and correlation (c) functional: 

771HSEO3 1 T7iHF,sr,u 1 ^ j7iPBE,sr,u , 77PBE,lr,u , 771PBE /-1 \ 

Acc = ^Ar + jAz + h x + ^0 > { l > 

where (sr) and (lr) refer to the short- and long-range parts of the respective exchange 
interactions, whereas fj, controls the range separation, varying between 0.2 and 0.3 A -1 . We 
have used \i = 0.2 A -1 . 

The Pd(100) supported c(4x2)-Mn 3 04 nanolayer has been modeled with a repeated 
slab constructed by a four layers thick Pd(100) substrate and a single manganese deficient 
MnO(lOO) overlayer with stoichiometry Mn 3 04. In order to allow for a rhombic distribution 
of manganese vacancies a two-dimensional (2D) c(4x2) unit cell has been utilised, which 
also permitted to test different magnetic configurations. Finally, particular care has been 
devoted to the interface stacking between the Mns04 overlayer and the metal substrate. 
Overall we have tested three models for the vertical set registry, which will be discussed in 
more details in the next session. 

As regards the structural relaxation, the two bottommost Pd layers were kept fixed. To 
relax the remaining atomic positions (the topmost 2 Pd layers plus the Mn 3 04 overlayer) 
we used the interatomic forces calculated through the Hellmann-Feynmann theorem and 
the geometry was optimised until the change in the total energy was smaller than 10~ 3 
eV between two consecutive ionic configurations. Well converged structural and electronic 
relaxations were reached at an energy cutoff of 300 eV and adopting a T-centred symmetry 
reduced 2D 6x6 Monkhorst-Pack [3J| fc-point mesh for Brillouin zone integrations. 

To determine the stability of the various models we have compared their generalised 
adsorption energy 7 defined as: 

7 — (-^slab — -Epd(lOO) — ^Mn^Mn ~ ^O^o) (2) 

where E s i a b indicates the total DFT energy of the Mn 3 O4/Pd(100) system considered, 
Epd(ioo) refers to the clean Pd(100) DFT energy, whereas n Mn / and /i Mn /o are the number 
of manganese (n=3) and oxygen (n=4) atoms and the corresponding reference energies. As 
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FIG. 1: (Colour online) Top view of the geometrical and magnetic models explored for the Pd(100) 
supported c(4x2) Mn304 phase. The models consist of a MnO(lOO) overlayer characterised by a 
rhombic distribution of Mn vacancies. Three different adsorption site registries of the Mn304 
layer on the Pd atoms underneath have been considered. Models RH1 and RH2 refer to a atop 
O-Pd and Mn-Pd adsorption configuration, respectively, whereas RH3 corresponds to a bridge site 
adsorption of both O and Mn on Pd(100). Plus and minus signs indicate the orientation of the 
Mn spins perpendicular to the surface, pointing inward and outward, respectively. Besides the FM 
arrangement two distinct AFM orderings were studied, which we have called AFM1 and AFM2. 



discussed elsewhere 



321 ] , we set f^un and fj>o to the energy of fee bulk Mn (7-Mn) and half of 



the free O2 dimer energy, respectively. 

Finally, phonon frequencies and eigenvectors have been evaluated by diagonalising the 
dynamical matrix of the most stable structure. 



III. RESULTS AND DISCUSSION 



The description of the structural and magnetic models considered in our study to simulate 
the experimentally observed c(4x2) phase are schematically depicted in FigflJ The common 
building blocks consist of a Pd(100) substrate containing 4 layers on which a compressed 
MnO(lOO) monolayer with a rhombic distribution of Mn vacancies, resulting in an Mn 3 C>4 
stoichiometry, has been placed. As already mentioned in the introduction the lateral com- 
pression compensates for the positive mismatch between MnO and Pd. We have therefore 
adopted for both substrate and interfacial layer the PBE minimised lattice constant ap^ E = 
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2.79. A, in good agreement with the experimental value of 2.75 A. Based on the square-like 
termination of MnO(lOO) and Pd(100) three models of the interface emerge naturally de- 
pending on the site registry: RH1 is constructed by placing the interfacial O on top of the 
Pd atoms; RH2 correspond to a top site Mn-Pd junction, whereas in the RH3 model both 
Mn and O are placed in the Pd(100) bridge sites. In terms of the Mn vacancies the above 
models correspond to a hollow-site (RH1), top-site (RH2) and bridge-site (RH3) vacancy. 
To investigate the effect of magnetism on the stability of the c(4x2)-Mn 3 O4/Pd(100) system 
we have computed the total energy of each structural model as a function of the Mn spins 
orientation. Besides the trivial ferromagnetic (FM) ordering, the choice of a c(4x2) unit 
cell allows for the set up of two additional antiferromagnetic (AFM) orientations, AFM1 
and AFM2. As displayed in FigJT], AFM1 and AFM2 differ for the dissimilar spin alignment 
along the [100] direction which is FM and AFM, respectively. 

Table [I] summarises the calculated relative energies for all nine possible configurations. 
The PBE values show that the RH1 is by 100-200 meV more favourable than any other 
geometrical setup and that the AFM2 ordering results in the most stable magnetic order. 
However, the inclusion of a fraction of non-local Hartee-Fock exchange within the HSE ap- 
proach, though favouring the RH1 geometry like PBE, reverts the magnetic relative stability 
preferring the FM configuration. Although the contradiction between PBE and HSE results 
is not surprising, in particular if related to magnetic properties of transition metal oxides, 
the small energy differences (tens of meV) and the lack of experimental magnetic analysis on 
this system do not permit a definitive, beyond doubt, conclusion regarding the most stable 
magnetic ordering. In order to arrive at an indisputable and trustworthy answer a more 
detailed investigation of both RH1-FM and RH1-AM2 models is required, which discloses 
the interplay between the magnetic ordering and the structural, electronic and vibrational 
properties for which accurate experimental findings exist. In the following we will therefore 
present in parallel the predicted physical properties of the FM and AFM2 RH1 models and 
compare them with available experimental results, with the aim of determining the most 
favourable magnetic ground state. 

In Fig. [2] we provide a graphical illustration of the optimised geometries for the RH1 FM 
and AFM2 models, whereas in Table [Til we list the most relevant structural relaxations and 
magnetic moments. Although the two minimised structures appear very similar, there are 
important magnetically driven differences which play a significant role, as we will discuss 
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TABLE I: Calculated relative stability (meV/unit cell) for the magnetic models investigated. HSE 
calculations have been only performed for the more favourable PBE models, namely RH1-FM and 
RH1-AFM2. 



FM AFM1 AFM2 

PBE 

RH1 41 97 

RH2 280 273 171 

RH3 230 301 180 

HSE 

RH1 - 22 



later on. First we focus our attention on the similar structural character of the explored 
models in terms of surface layer displacement (8), vertical buckling (b) and interface Mn 3 04- 
Pd(100) distance (z). Most of the relaxations take place in the Mn 3 04 layer which is found 
to be well separated from the substrate by z?»2.3 A, about 20% larger than the Pd bulk 
interlayer distance. As a consequence, the Pd(100) substrate remains structurally unaffected 
by the growth of the Mns04 monolayer. The presence of the Mn vacancies dominates the 
structural rearrangement and determines the appearance of two distinct Mn types: Mnl 
localised between two nearest neighbour vacancies and Mn2 forming zig-zag Mn-Mn chains 
along the [100] direction. In addition, the formation of the Mn vacancy breaks the local Mn- 
O bonds and induces a planar outward relaxation of the oxygen atoms which move closer 
to Mnl, resulting in the waving long-range arrangements highlighted in Fig. [2j The lateral 
oxygen displacement 5o=0.12 A is found to be identical for both FM and AFM2, whereas 
the movements of the manganese atoms is notably different. Within an AFM2 ordering 
the lateral Mn2 displacement (<5mii2 = 0.12 A) turns out to be 0.04 A smaller than in FM 
and as a consequence of that a small corrugation of 0.05 A along the chain is observed in 
RH1-AFM2 between adjacent antiferromagnetically coupled Mn2 atoms. Finally the Mnl 
atoms in RH1-AFM2 experience a shift (#Mni=0.07A) towards the nearest AFM coupled 
Mn2, which is responsible for the smaller width of the Mn2-Mn2 zig-zag chain in RH1- 
AFM2 and the quite large buckling of 0.33 A between inequivalent oxygen atoms. Within 
the FM configuration the 8mui movement is forbidden by symmetry and cannot take place, 
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FIG. 2: (Colour online) Optimised geometrical structure for the c(4x2) Mn304 RH1 model in FM 
(left) and AFM2 (right) magnetic configuration. The Pd(100) substrate is displayed with light 
gray large spheres, whereas Mn and O atoms are depicted with small dark (blue) and gray (pink) 
spheres. Two distinct Mn species are distinguishable: Mnl, sandwiched between two vacancies 
and two Mn2 atoms forming the zig-zag arrangements highlighted by the horizontal lines. Atoms 
around the manganese vacancies experience considerable strains (indicated by arrows) which leads 
to similarly distorted FM and AFM2 structures. The only sizeable differences between the FM and 
AFM2 minimised geometries reside in the displacement <5mh1 an d in the vertical buckling between 
oxygen atoms which are not present in the FM model. 

thus preventing the corrugation between the oxygen species. 

As for the local magnetic moments (also reported in Table [TT]) we do not observe any 
significant difference between the two spin arrangements. The inclusion of a fraction of HF 
exchange increases the local Mn moments by 10-15% with respect to the PBE values and 
gives an induced magnetic moment of ~ 0.4 \xb on the Pd atoms directly underneath the O 
atoms. These HSE results are not unexpected considering the tendency of hybrid functionals 
to provide a generally more localised description and a larger splitting between majority and 
minority states, as shown in the comparison between PBE and HSE Mn density of states 
(DOS) in Fig. El 
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TABLE II: Calculated optimised geometry and local magnetic moments of the Mn30 4 RH1 phase 
compared with the ideal Mn404 structure, for both FM and AFM2 magnetic orderings. The 
structural data refer to PBE calculations and are all given in A. Below the PBE magnetic moments, 
in fis, the corresponding HSE values are also listed. Symbols (schematically explained in Fig. [2]): 
5 refers to the in-plane displacements of Mn and O atoms; b's are the vertical bucklings and z 
indicates the average interlayer distance between the MnO layer and the Pd(100) substrate (also 
given in percent relative to the bulk Pd(100) interlayer distance). 
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To understand the stabilisation mechanism of the RH1 phase it is instructive to make 
a comparison with the stoichiometric MnO monolayer. The epitaxial growth of a perfect 
MnO(100) monolayer on Pd(100) is clearly unfavoured because of the large lattice mismatch 
between MnO and Pd. In fact, the absence of Mn vacancies prevents the planar relaxations 
which give rise to the RH1 structure and lead to a very corrugated layer with a buckling of 
0.5 A between O and Mn as reported in Table ILT1 In addition, for the ideal MnO(100) system 
we predict a larger separation z between adlayer and substrate which weakens substantially 
their mutual interaction. 

Coming back to the nonstoichiometric adlayer, the analysis of the DOS combined with 
the structural properties described above allow for a comprehensive understanding of the 
bonding picture. The hybridisation between Pd d z i and O p z orbitals, particularly strong 
around the Fermi energy (Ep), reflects the direct Pd-0 atop site adsorption specific of the 
RH1 model. The electronic interactions within the Mn 3 4 layer are instead dominated 
by the planar bonds between Mnl and Mn2 d xy states with O p x and p y orbitals. The 



9 



FM 



AFM2 




An 




Mn1,d 




! — v 

— 








/y 




1 i 

4h — 1 — ^ 




<\/r 


i i i 
Mn2, d 


— i — i — i — 

A. , 




~w — ' 


1 


711 




f 


-M — -+- 






i i i 
0,p 


— i — i — i — 


Pz " 








hr 








1 — 1 — (— 

' ^surf, d 


— i — i i , 


— i — i — i- 

d z ,- 






r " 







0.5 



0.3 


-0.3 
0.2 


-0.2 

0.3 



-0.3 



-6 -4 -2 

Energy 



2 4 6 -6 -4 -2 2 4 6 

[eV) Energy (eV) 



-6 -4-2 2 4 6 -6-4 -2 2 4 6 

Energy (eV) Energy (eV) 



FIG. 3: (Colour online) Atom and spin resolved local density of states (LDOS) of the RH1 structure 
calculated within the HSE approach in the FM and AFM2 magnetic configurations. Shadow areas 
(left side of panels) refer to PBE data. For each atom, majority and minority spin states are 
plotted in different panels, as indicated by arrows. 



structural differences between Mnl and Mn2 are clearly expressed by a dissimilar DOS. More 
specifically, the Mnl-0 bond, broken by the creation of the Mn vacancy, yields the formation 
of well localised states at Ep in Mnl which are not observed in Mn2. These localised d 
states which characterise both FM and AFM2 magnetic alignments are responsible for the 
larger spin splitting and magnetic moment in Mnl (mMni ~ 4.2 ixb) as compared to Mn2 
(^Mn2 ~ 3.7 (j,b)- A second important difference between FM and AFM2 resides in the 
DOS near Ep, in the energy window between -2 and +2 eV. In the FM alignment this 
energy region shows a predominant Mnl d xy spin up nature, characterised by a 3 eV broad 
band crossing Ep, and only above 1 eV the Mn2 contribution becomes significant. Within 
the AFM2 ordering the Mnl d xy band near Ep is much more localised and displays a single 
peak structure which decays sharply at Ep. Since the Mn2 d states remain unchanged as 
compared to the FM situation, a gap region emerges in the range (0-0.7) eV where the 
Mns04 overlayer shows an electronic behaviour different from that predicted for the FM 
phase. 

At this point, comparison with the experiments is required to validate our results and 
to finally address the issue of which of the two modelled magnetic structures corresponds 
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FIG. 4: (Colour online) HSE RH1-FM and RH1-AFM2 simulated STM images compared with the 
experimental STM image (shown in the insets) taken at a sample bias U=+0.8 V and tunnelling 
current I = 0.25 nA. Oxygen atoms are not "seen", whereas manganese atoms appear as bright 
maxima. The detailed analysis of the maxima is discussed in the text. 

to the observed one. In FigJH we compare the calculated FM and AFM2 STM images 
with the experimental picture. Remarkably, the comparison unambiguously shows that the 
theoretical RH1-FM model resembles much better the experimental STM image as far as the 
contrast and the relative brightness of the observed maxima are concerned. The measured 
picture is characterised by a rhombic arrangement of bright spots and dark depressions 
connected by weaker segments. These features are perfectly well reproduced by the simulated 
FM STM picture and can be interpreted in the following way: the black depressions reflect 
the network of Mn vacancies, the single bright spots can be assigned to the Mnl species 
whereas the light segments correspond to the zig-zagging Mn2-Mn2 chains embedded in the 
regular array of oxygen atoms, which are not seen in the STM image. In the AFM2 model 
the Mnl and Mn2 relative brightness is reverted: The zig-zagging Mn2-Mn2 chains appear as 
the prominent feature, whereas the Mnl species are seen as weaker unresolved connections. 
Now, the following question arises naturally: why do the FM and AFM2 simulated profiles 
display different features? To answer this issue we note that the experimental STM image 
has been taken at bias +0.8 V, for which the FM and AFM2 phases show a very dissimilar 
character as discussed above and emphasised by the vertical dashed lines in Fig. [3j In 
the FM structure the STM simulation reproduces the states on top of the Mnl d xy band 
centred around Ep weakly overlapping with Mn2-like states. In contrast, within the AFM2 
configuration the only states contributing to the DOS near +0.8 eV come from the Mn2 d xy 
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TABLE III: Predicted dipole active modes for the RH1 phase within the FM and AFM2 magnetic 
orderings compared to the measured HREELS phonon peak. All data are given in meV. 





PBE 


HSE 


Expt. 


FM 


36.4 


43.3 


43.5 


AFM2 


35.4 


39.6 





orbitals, which are therefore the ones seen in the AFM2 STM, in line with the interpretation 
given above. 

Summing up, the distinct electronic character of the FM and AFM2 phases near the 
experimental bias voltage is correctly recognised by HSE and is reflected in the simulated 
STM images. The FM alignment of spins is essential to reproduce the experimental STM 
features and the HSE approach, by reverting the PBE relative stability between the FM 
and AFM2 structures (see Tabled]), shows once more its predictive power and its capabil- 
ity of improving the description of transition metal compounds with respect to standard 
DFT. We emphasise that a standard PBE approach alone cannot provide a satisfactory and 
comprehensive understanding of this system. 

Additional support for the proposed RH1-FM surface structure is provided by the com- 
parison between experimental and theoretical vibrational properties reported in Table IHIi 
HREELS experiments established that the phonon loss spectrum is characterised by a sin- 
gle peak located at 43.5 mev[l^]. Indeed, our calculations correctly predict a single active 
dipole mode and demonstrate that the FM structure perfectly reproduces the experimen- 
tal vibrational energy of 43.5 meV. Two important conclusions deserve particular enphasis: 
First, PBE is unable to supply a satisfactory answer, providing a much too low dipole active 
mode (~ 36 meV) and second, the AFM2 phonon peak is found to be 4 meV lower than 
the measured value. To conclude, in Figj5]we show a pictorial view of the atomic displace- 
ments associated with the vibrational peak. Upon diagonalisation of the dynamical matrix 
HSE establishes that the active dipole mode originates from the anti-phase and out-of-plane 
vibrations of O and Mn atoms, a vibrational mode typical of the ideal MnO(lOO) surface loj ]. 

These results as well as the explanation of the STM findings are the most convincing 
evidence for the correctness of our computational model. 
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FIG. 5: Schematic illustration of the atomic displacements giving rise to the calculated Mn-0 
vibration mode. Mn and O atoms are sketched as black (blue) and gray (pink) spheres, respectively. 

IV. CONCLUSIONS 

The formation of the two-dimensional c(4x2)-Mn 3 04 nanolayer phase on Pd(100) has 
been explored by means of first principles PBE and HSE calculations. Combining struc- 
tural, electronic, magnetic and vibrational analysis we have demonstrated that the preferred 
model which correctly describes the available experimental information entails the formation 
of rhombically distributed Mn vacancies in a compressed MnO(lOO) adlayer, preferentially 
adsorbed with an O-Pd atop site registry on the Pd(100) substrate. We found that hybrid 
functional methods such as HSE are essential to capture and explain all relevant proper- 
ties, which would be otherwise only loosely understood within a standard DFT approach. In 
particular, HSE reverts the relative stability of the two most stable magnetic structures con- 
sidered and improves the computed vibrational energies over PBE. The role of the magnetic 
interactions is crucial to interpret the STM and HREELS experiments. The AFM2 ordering 
describes incorrectly the DOS of the Mns04 adlayer, reverts the brightness of the observed 
STM maxima and decreases significantly the phonon dipole active mode, in contrast with 
the FM results which are fully consistent with the experimental observations. 

Our study underlines the complex interplay between orbital, spin and lattice degrees of 
freedom occurring in 2D transition metal oxide nanolayers and shows that the post-DFT 
approach we have used represents a qualitative and quantitative step forward in surface 
science modelling. 
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